Holography and the speed of sound at high temperatures 



Paul M. HohleiB and Mikhail A. Stephano\jl| 

Department of Physics, University of Illinois, Chicago, IL 60607-7059, USA 

We show that in a general class of strongly interacting theories at high temperatures the speed of 
sound approaches the conformal value = 1/3 universally from below. This class includes theories 
holographically dual to a theory of gravity coupled to a single scalar field, representing the operator 
of the scale anomaly. 
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Introduction — The discovery of the holographic cor- 
respondence between string and gauge theories [H, 0, S| 
has led to a number of applications to modelling proper- 
ties of quark-gluon plasma of QCD in the regime relevant 
to heavy-ion collision experiments (see, e.g., Ref. Q for 
a guide and references). 

Although the precise form of the theory holographi- 
cally dual to QCD is not known, many features of QCD 
can be implemented on the string theory side of the cor- 
respondence. Alternatively, one can pursue a bottom-up 
approach, by constructing an effective theory, or a model, 
with a minimal set of operators needed to describe the 
relevant physics. In this paper, we consider a general 
class of such models, describing the physics of scaling vi- 
olation in QCD. In Refs. [1, @, [3| such theories have been 
considered with the aim to model, or mimic, the equation 
of state of QCD known numerically from the lattice cal- 
culations. This has been achieved by tuning the potential 
of the scalar field in the model. 

The main goal of this paper is to establish properties 
of the equation of state p(e), which are universal with 
respect to the choice of the potential V{(j)). An obvi- 
ous universal property is that the conformal symmetry 
is gradually restored, and thus p{e) — *■ e/3, as the tem- 
perature T oo. This is natural since the temperature 
becomes the only relevant scale. It is a welcomed fea- 
ture of the models since the same phenomenon occurs in 
QCD. 

Here, we shall show that the deviation of the speed of 
sound = dp/de from the conformal limit = 1/3 is 
universally negative in such holographic models, at least 
to the leading order in inverse temperature. 

At the outset, we must point out that — 1/3 is by 
no means a universal upper bound on the sound veloc- 
ity. A counterexample is given by the speed of sound 
in QCD at large isospin chemical potential [sl, where Cs 
can approach the speed of light in a certain limit. Never- 
theless, a number of string theory examples ^, 1^, 11 1 of 
holographically dual theories do indeed consistently show 

< 1/3. 

In QCD, Monte Carlo lattice calculations [l2] show 
that this inequality is fulfilled. In the regime T 3> Aqcd, 
this inequality follows from asymptotic freedom: /3(a) < 
0. For a pure A^c-color Yang-Mills theory at high T, for 
example, 13] 



How general is the inequality < 1/3 and what are 
the prerequisites for it to hold? For example, does it hold 
in regimes where the theory is close to being conformal, 
but is still strongly interacting, as is apparently the case 
for QCD in the range of T - (2 - 3) TJ 

To investigate this question, we consider a class of the- 
ories, or models, which possess a holographically dual 
description. The minimal set of operators that we need 
to consider consists of the stress-energy tensor T^^"^ and 
a scalar operator O necessary to provide a nonzero right- 
hand side of the scale (trace) anomaly equation. As an 
example, in QCD, the role of such an operator is played 
by TrF^: = P{a)/{8Tra^)TiF^. 

At this point it is helpful to express the energy density 
e and pressure p in terms of the heat function (enthalpy) 
w = e + p and the trace of the stress-energy tensor = 
e — 3p: 



Then 



e = (3w + 9)/4; p ^ {w 6')/4. 



dp 1 — dO/dw 



de 3 -f dO/dw 



(2) 



(3) 



In conformal theories 9 = Q and = 1/3. The inequal- 
ity < 1/3 is equivalent to dO/dvu > 0. This inequality 
is somewhat reminiscent, but not equivalent, to the in- 
equality d{pT-'^)/dT = 0/T^ >0 conjectured in Ref.[li| 
and the weaker constraint JdT6/T^ > proposed in 
Refs. [ll,[li. 

Since the enthalpy w appears to play here the role of 
a more natural thermal parameter than T, it is helpful 
to keep in mind that w is a monotonous function of T: 
dw/dT = Cy + s > 0. 

The theory — The class of the five-dimensional (5D) 
holographic theories that we consider contains a scalar 
field (j), holographically dual to the (scaling violating) op- 
erator O, coupled to the metric (dual to T^'^). The 
fields are governed by the action 




(4) 



~ 1/3 + 5A^c/(367r) /3(a) < 1/3. 



(1) 



where R is the Ricci scalar, g is the determinant of the 
metric, 7 is the determinant of the induced metric on the 
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UV boundary 9M, K is the extrinsic curvature on 9M, 
and is the 5D Einstein gravitational constant. The 
value of is inversely proportional to the number of the 
degrees of freedom in the dual four-dimensional theory, 
e.g., iV^ in a gauge theory with large number of colors Nc- 
The smallness of (i.e., the largeness of the number of 
colors) controls the semiclassical approximation which we 
use. The last term in the action is the Gibbons-Hawking 
term, which removes the boundary terms arising upon 
integration by parts of the terms in R linear in second 
derivatives of the metric [l7| . This boundary term does 
not affect classical equations of motion, but is essential 
for evaluating variations of the action with respect to the 
boundary values of the fields. 

We assume that the potential V{4i) has an extremum 
at = 0. The value of ^(0) = 2A = -12, where A 
is the cosmological constant necessary to achieve proper 
asymptotics of the metric near the boundary. 

It should be noted that we are going to treat the holo- 
graphic theory with a method technically different from 
the one used in either Refs. [1, @] or Ref. Our ap- 
proach is complementary to the existing ones and has 
its own advantages, which might go beyond the specific 
application which we consider in this paper. 

According to the holographic correspondence, the cor- 
relation functions of the dual four-dimensional (4D) the- 
ory are equal to the variations of the 5D action under 
the changes of the boundary conditions on the 5D fields. 
We shall first determine the extremum of the action, and 
then consider first-order variations, which we can relate 
to energy density and pressure. The most general met- 
ric (up to general coordinate transformations) possessing 
three-dimensional (3D) Euclidean isometry is given by 

ds^^\{-f{z)dt' + dx^)+e^''^'-^^^^. (5) 

The functions B{z) and f{z) depend on the holographic 
coordinate z and will be determined by extremizing the 
action. The boundary dM is located at z = e with e 
acting as an UV regulator. In relations involving only 
physical quantities e can be taken to 0. 

Substituting the metric of Eq. ([5]) into the equations 
of motion wc find: 

B = -l^\ (6) 

/ = (4 + s) ./, (7) 

-6/ + / (24 - </)2) + 2e2-s 1/(0) = 0, (8) 

^/ + 0(/-/(4 + i?)) -e2^r(0) = O, (9) 

where a dot denotes a logz derivative, e.g., = zd(j)/dz, 
while V' = dV/dcf). The holographic correspondence pro- 
vides the boundary conditions at the UV boundary z ~ e. 
The boundary condition on 0, 

0(e) =c£^-, (10) 



corresponds to introducing a term c O into the action of 
the 4D theory dual to the theory in Eq. (|4|), where c is 
the source of the scalar operator O, dual to the field 0. 
The dimensions of the source and the operator are given 
by 

[c] = A_; [0]=A+; A+ + A_ = 4. (11) 

The metric must approach the Lorentz invariant form 
at the UV boundary, hence, 

fie) = l- (12) 

Equation ([7]) can be integrated once to give 

/ = -wz^e^. (13) 

The integration constant w must be positive if the metric 
is to possess a horizon /{zh) = at some value of zh- 

There are no more independent boundary conditions. 
The boundary value of B is determined by Eq. ([8|) , which 
is algebraic in B. The role of the second boundary con- 
dition for Eq. ^ is played by the requirement that is 
finite at the horizon, z = zh, which is a regular singular 
point of the second order differential Eq. ([9]) . 

Thus, we find a two-parameter family of solutions. 
These parameters are c from Eq. pO|l and w from 
Eq. (fT3|) . Different members of this family are related 
by rescaling z z/X, w X'^w, c — > A^^c, which rep- 
resents scale invariance inherent in the action given by 
Eq. Q. 

Near the z = e — > boundary, 0— >0, B— >0, -B— >0 
and / ^ 0. Equation ([9|) for can be linearized and 
the asymptotic behavior of near the boundary can be 
determined easily: 

0(z) (c - de^+-^-)z^- + dz^+ -f . . . , (14) 

where the curvature of the potential V"{0) = deter- 
mines the indices A± = 2 ± VT+rr?. The coefficient of 
the first term is related to c by Eq. (fTO|) . The coefficient 
d of the second linearly independent solution should be 
determined by the finiteness condition at the horizon. 

By calculating the derivative of the 5D action with 
respect to c and matching it, by holographic correspon- 
dence, to the expectation value (O), one finds 18] 

(O) = -9S'5/9c= -d(A+ - A_). (15) 

Equation of state — Rather than following existing ap- 
proaches to calculating the speed of sound, based on the 
relation = d log T/d logs, or looking at the poles 
of the two-point function for the sound channel [HI [2^ , 
we shall directly calculate owe-point functions (T**) and 
(T^^) and use the relation = dp/de. As we shall see, 
this quickly yields some closed-form results, not evident 
in the alternative approaches. 

In order to calculate the one-point function, such as 
(T**), we observe that the generalization of the 4D theory 



3 



to a curved 4D background with metric /i^^ corresponds, 
holographically, to imposing the following boundary con- 
dition on the bulk 5D metric 



Holographic correspondence then gives us 



2 (555/(5 V- 



(16) 



(17) 



Variation of the boundary condition causes variations of 
the metric in the bulk, as it follows the equations of mo- 
tion. But since we are varying around the cxtrcmum 
of the action, the only contribution to the first-order 
variation comes from the boundary. This contribution 
consists of the variation of the boundary terms, which 
appear during integration by parts while deriving equa- 
tions of motion for, e.g., gtt, as well as the variation of 
the Gibbons-Hawking term. Altogether this gives 



w - 



(18) 



Both of these quantities are divergent as e 0. This 
represents the familiar vacuum energy divergence in the 
quantum field theory. A simple vacuum subtraction 



(r")-(r")T=o, 



(19) 



T=0, 



takes care of this, and we find for the energy density and 
the pressure, after solving Eq. (I8| for J5 at z = £ — + 
with (j) from Eq. (fT4|) : 



(20) 



e = (3w-cdA_ (A+ - A_))/4, 
p = (w + cdA- (A+ - A_))/4. 



We can now conclude that the integration constant in- 
troduced in Eq. (fT5|) is the enthalpy, w = e + p, of the 
corresponding 4D field theory. We also observe that the 
scale anomaly 



3p: 



-cdA_(A+- A_) = A_c(0), (21) 



is related to the expectation value of the operator O [us- 
ing Eq. (Uni)]. This relation is easy to derive directly on 
the field theory side by observing that the violation of 
the scale invariance comes from the dimensionful param- 
eter c. Naturally, it is proportional to the dimension of 
c, A_. 

From Eq. (^0]) . the speed of sound can be expressed 
now in terms of the derivative dd/dw at fixed c and di- 
mensions A-I-: 



dp _ 1 -|-cA_ (A+ - A_) {dd/dw) 
d^ ~ 3 - c A_ (A+ - A_) [dd/dw] 



(22) 



It should be emphasized that this equation is an exact 
expression valid for all temperatures. 



High temperature — The high-temperature limit is 
equivalent to the small c limit, since c is the only 
other dimensionful parameter in the theory, and high- 
temperature expansion is controlled by the dimensionless 
parameter c/T^", or cw~^~/'^. In this limit. 



1/3 + (4/9) cA_ (A+ - A_) [dd/dw) + 



(23) 



We shall now show that, even though generally the de- 
pendence of d on w can be only found numerically, in the 
high-temperature (i.e., high w) limit an analytic expres- 
sion can be found for arbitrary potential V[(j)). 

One can begin by observing that at large w the func- 
tion / varies very rapidly according to Eq. ()13p . This 
means one can neglect variations of the function B be- 
tween the boundary z = e and the horizon z — zh, since 
zh becomes small (as w~^^^). Since on the boundary 
B = (up to terms of order e^'^- , negligible here, ac- 
cording to Eq. dH)), we find from Eq. (fT3|) 



f[z) 1 -wz'^/A. 



(24) 



Another consequence is that (j), which is small at z = s, 
remains small up to zh [(f> ~ c/T^- <^ 1), and 

the linearized approximation to Eq. ([9|) is valid not only 
near the boundary, but all the way to the horizon. With 
B = and / from Eq. (|24|) we obtain 



1 — - WZ' 



3 WZ 
-z + — 



= 0. (25) 



Equation ([25|1 can be solved analytically 

0(z) = cz^- 2F1 (A_/4, A_/4, A_/2, wz^ A) 
+ dz^+ 2F1 (A+/4, A+/4, A+/2, wz^A) , 



(26) 



where the coefficients follow the notations of Eq. ([T 
(up to terms 0[e^+~'^-), here negligible). Both lin- 
early independent solutions are logarithmically divergent 
at the horizon z = zh, where wzj^/A = 1. The condition 
|(/)(z^f)| < 00 requires us to select the linear combination 
in which these divergences cancel. This fixes d in terms 
of c: 

d = -cw;(^+-^-)/4i?(A_), (27) 
where the function £'(A_) = 1/£'(A+) is given by 

D[A^) 



'^'^ cot(.A_/4)S^^. (28) 



A_ 



r(A_/4)4 



Substituting Eq. (|27p into Eq. we find our main re- 
sult 



c^ = i-ic^A_(A+-A_)^«.-^-/^i?(A_ 



(29) 



It is clear that the correction term is negative for all val- 
ues of < A_ < 2 (i.e., 2 < A+ < 4). As expected, the 
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correction vanishes with the power of w dictated by the 
dimension of the operator O. 

Discussion — When this work was completed, the au- 
thors learned about a similar result [2l|, obtained using 
the relation = d log T/d logs. In Ref. [2l| the speed of 
sound is expressed in terms of the value (j)H of the scalar 
field at the horizon (since both s and T are calculated at 
the horizon): = 1/3— C(j>j^ + . . . . The coefficient C is 
then given by an integral from the boundary to the hori- 
zon of the square of a hypergeometric function, which is 
manifestly positive. We have verified that, upon evalu- 
ating the integral, the result of Ref. 2l| coincides with 
Eq. ([29|) . For completeness of comparison we evaluate 
explicitly by using Eq. (|26p and applying relation 
Eq. dSID 



OH 



c w 



A/4-1 2-^/2(2A-4) 



r(A/4)2 



1 

18^ 



(30) 



(4- A) (4-2A) tan (7rA/4) 4>jj + . . . . (31) 



where A — A^. Using Eq. ((29|) we then find 



As we have pointed out already, and as Eq. ([50]) 
demonstrates explicitly, in the high-temperature 
(large w) limit the value of remains small everywhere 
between the boundary and the horizon. This is the 
origin of the universality we find. The speed of sound 
near the high T limit depends only on the behavior of 
the scalar potential near the origin, i.e., specifically, on 



V"{0) 



A(A-4). 



One can interpret and further generalize our results 
in the following way. In the high-temperature limit, 
Eq. (PT|) becomes 9 — A_ xo + ■ ■ ■ ■ The susceptibility 
Xo = d{0) /dc can be related to the effective potential 
for (O), defined as the Legendre transform of the gen- 
erating functional W{c) (in holography W = —S5), i.e., 
r((0)) = W{c) + c{0), as xo = l/r"(0). stability im- 
plies r"(0) > 0, and thus xo > 0. Consequently 6* > 0, 
as conjectured in [l3 | . Holographic models confirm these 
expectations according to Eqs. (fT5|) and (p7| . 

We also find that dO/dw > 0, which requires 
dxo/dw > 0. That means that the curvature r"(0) = 
Xo^ decreases with temperature. This behavior is un- 
usual, if one recalls that in a weakly coupled scalar 
theory, the leading perturbative temperature correction 
to the curvature (~ XT^) increases with temperature 
(e.g., the restoration of a broken symmetry is a well- 
known manifestation of this). The opposite behavior of 
the curvature of the effective potential r{{0)) can be un- 
derstood by counting dimensions: [xq^] — A_ — A+ < 0, 
which means Xo^ ~ j-a^-a^ ^^^^ decreases with T. 

In conclusion, we have shown that in a quite general 
class of gravity dual theories with a single scalar operator 
representing the scale anomaly the speed of sound always 
approaches the conformal value = 1/3 from below. 
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